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A REEVALUATION OF DATA ON ATMOSPHERIC TURBULENCE AND AIRPLANE GUST LOADS 

FOR APPLICATION IN SPECTRAL CALCULATIONS > 

By Harry Press, May T. Meadows, and Ivan Hadlock 


SUMMARY 

The available information on the spectrum of atmospheric 
turbulence is first briefly reviewed. On the basis of these results, 
methods are developed for the conversion of available gust statis- 
tics normally given in terms of counts of gusts or acceleration 
peaks into a form appropriate for use in spectral calculations . 
The fundamental quantity for this purpose appears to be the 
probability distribution of the root-mean-square gust velocity. 
Estimates of this distribution are derived from data for a number 
of load histories of transport operations; also, estimates of the 
variation of this distribution with altitude and weather condition 
are derived from available data and the method of applying these 
results to the calculation of airplane gust-response histories in 
operations is also outlined. 

INTRODUCTION 

During the last few years, advances have been made in the 
analysis of airplane behavior in rough air through the appli- 
cation of the theory of random processes and the techniques 
of generalized harmonic analysis (refs. 1 to 7) . The applica- 
tion of these techniques is based upon the representation of 
atmospheric turbulence as a continuous random disturbance 
characterized by power-spectral-density functions and certain 
probability distributions. The power spectrum of the tur- 
bulence is then used along with the airplane response char- 
acteristics to determine the power spectra and other statis- 
tical characteristics of the airplane loads or motions in rough 
air. The application of this approach to the problem of 
calculating load and other histories for operational flight 
requires detailed information on the spectrum of turbulence 
in the atmosphere. The information required may be con- 
sidered of two types: detailed information on the spectrum 
of turbulence and its variations, and information on the 
probability of encountering the various spectra. 

A few measurements of the power spectrum of atmospheric 
turbulence have so far been made and are briefly reviewed in 
the report. These measurements indicate that, over most of 
the frequency range of interest, the spectra are inversely 
proportional to the square of the frequency and may be 
approximated by simple analytic expressions such as have 
been used in wind-tunnel studies of isotropic turbulence. 


Also, the intensity of the turbulence as given by the root- 
mean-square gust velocity varied appreciably with the 
weather conditions. These results thus appear to provide 
some information on the power spectrum and its variations. 
They do not, however, provide any information of the second 
type required, that is, information on the probability of en- 
countering the various spectra in actual operations. The 
purpose of this report is to provide information of this second 
type, that is, information on the probabilities of encountering 
the various conditions of turbulence. 

The only source of information on the probability of en- 
countering the various conditions of atmospheric turbulence 
appears to be the considerable amount of statistical data con- 
cerning atmospheric turbulence and airplane loads in rough 
air that has been collected by the National Advisory Com- 
mittee for Aeronautics in the last 20 years. (See, for example, 
refs. 8 to 12.) These data have been obtained, in most 
cases, from airplane acceleration measurements in normal 
operations, although in some instances (ref. 10) the data were 
obtained in special flight investigations. These data are 
generally given in the form of the number of peak accelera- 
tions or effective (or “derived”) gust velocities per second 
which exceed given values and in this form do not appear 
applicable to spectral methods of analysis. 

Fortunately, in the theory of random processes, relations 
have been derived (ref. 13) between peak counts (such as 
have been made for normal acceleration) and the associated 
power spectra. These relations apply to the case of a 
stationary Gaussian random process, the stationarity of 
the process implying that its characteristics do not change 
with time and the term Gaussian designating a process 
characterized by a Gaussian probability distribution for 
the amplitude of the disturbance as well as for its time 
derivatives. The approximately Gaussian character of 
turbulent velocity fluctuations has been noted and, for the 
case of atmospheric turbulence, results reported, for example 
in reference 3, appear to support such an assumption. 
Inasmuch as the intensity of turbulence is known to vary 
widely with weather conditions, the overall gust and load 
experience in operations cannot be considered a stationary 
Gaussian process. If the operational gust history is 


1 Supersedes NACA Technical Notes 3362 by Harry Press, May T. Meadows, and Ivan Hadloek, 1955, and 3540 by Harry Press and May T. Meadows, 1955. 
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considered to be a nonstationary Gaussian process varying 
only in intensity or root-mean-square gust velocity, the 
problem of specifying the gust history is reduced to that 
of specifying the probability distribution of the root-mean- 
square gust velocity. For this case, it appears possible to 
extend the results of reference 13 for the stationary Gaussian 
process in order to derive a basis for estimating the distri- 
bution of root-mean-square gust velocity from counts of 
acceleration peaks. 

On the basis of the foregoing considerations, techniques 
are developed in this report for the estimation of the prob- 
ability distributions of the root-mean-square acceleration 
and root-mean-square gust velocity from data on peak 
gust accelerations. Two approaches are considered: first, 
that the root-mean-square gust velocity takes on only a 
few discrete values and, secondly and more realistically, 
that the root-mean-square values cover a continuous range 
of values. These techniques are them applied to available 
gust-load data obtained from a number of transport opera- 
tions and th^ associated probability distributions of root- 
mean-square gust velocity are derived. Also, available 
data are used to estimate the variations in this distribution 
with altitude and weather condition. The method of 
application of these results to the calculation of gust-load 
and other airplane response histories in operations is out- 
lined and some of the limitations of the present results are 
indicated. 
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SPECTRAL CHARACTERISTICS OF ATMOSPHERIC 
TURBULENCE 
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In the way of background for the present analysis, it will 
be helpful to review the available measurements of the 
power spectrum of atmospheric turbulence. Figure 1 con- 
tains a summary of most of the available airplane measure- 
ments of the power spectrum of atmospheric turbulence. 
The first of these measurements was made by Clementson 
at the Massachusetts Institute of Technology (ref. 4), and 
subsequent measurements were made by the National Ad- 
visory Committee for Aeronautics (refs. 5 and 7), Douglas 
Aircraft Company, Inc. (ref. 14), Massachusetts Institute 
of Technology (ref. 15), and Cornell Aeronautical Laboratory 
(refs. 16 and 17). The curves shown represent the various 
power spectra which were obtained under different weather 
conditions. The abscissa is the frequency argument Q, which 
has the dimensions of radians per foot and is equal to 2?r 
divided by X, the gust wavelength. (The data shown cover 
a range of gust wavelengths from about 10 feet to 5,000 feet.) 

The spectra in all but one case are for the vertical or lateral 
component of the turbulence. In one case marked by the 
letter H, the spectrum is for the horizontal or longitudinal 
component of turbulence. Examination of these results indi- 
cates that the spectral shapes appear to be relatively consist- 
ent; in all cases, the power decreases rapidly with increasing 
frequency. In fact, in most cases, the spectra appear to be 
inversely proportional to the square of the frequency. This 
spectral shape of I/O 2 is in reasonable agreement with 
theoretical results obtained for the spectral shape at the 
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Figure 1. — Summary of airplane measurements of the power speetrum of atmospheric turbulence. 

higher frequencies in the theory of isotropic turbulence. At 
the lower frequencies, the situation is not as clear, few 
measurements being available for frequencies £2<C 0.002. 

Some additional measurements obtained at the Cornell 
Aeronautical Laboratory and at the NACA (ref. 18) do 
cover these lower frequencies and indicate a flattening of 
the spectrum at frequencies ft<C0.001. 

In addition to these variations in spectral shape, the 
various measurements also differ in turbulence intensity. 

The individual root-mean-square values are estimated to 
vary from roughly 1.5 to perhaps 8 feet per second, which, 
as will be seen, represent the relatively light-to-moderate 
levels of atmospheric turbulence. 

Another source of information on the spectral character- 
istics of atmospheric turbulence is the measurements at 
lower altitudes made from meteorological towers. A large 
number of such spectral measurements have now been 
obtained. A few representative measurements obtained at 
an elevation of about 300 feet and for various conditions of 
average wind speed V w (ref. 19) are shown in figure 2. 

These measurements extend to lower frequencies (longer gust 
wavelengths) than do most of the available airplane measure- 
ments. At the higher frequencies, these results approximate 
the same form 1/S2 2 that is characteristic of the airplane 



Figure 2. — Measurements of the power spectrum of atmospheric 
turbulence obtained at 300 feet from a meteorological tower (ref. 19). 
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measurements. In addition, at the lower frequencies, a 
definite tendency toward a flattening of the spectra can be 
noted. The variation in the spectrum intensity with wind 
speed V w should also be noted. 

Because of the general characteristics of these spectral 
measurements, it has *been convenient in theoretical studies 
to use the following analytical expression for the turbulence 
spectrum 




L 1+3 tt 2 L 2 

X (1 + S2 2 L 2 ) 2 


( 1 ) 


For design purposes, the overall gust experience in opera- 
tions is of concern. Presumably, the overall experience 
consists of various exposure times to each of the spectra 
shown in figures 1 and 2 as well as to other spectra associated 
with different weather conditions. It thus appears important 
to determine the proportion of flight time spent under these 
various conditions of turbulence. In the neit section 
methods are developed for obtaining information on this 
problem from available gust statistics. These methods are 
then applied in a subsequent section. 


where 12=2x/X and X is the gust wavelength. This expression 
has been useful in wind-tunnel studies of turbulence and 
has the general characteristics of the measured spectra of 
atmospheric turbulence. The equation has two parameters: 
the mean-square gust velocity ov 2 which describes the overall 
intensity, and the so-called scale of turbulence L which, in a 
sense, can be considered to be proportional to the average 
eddy size. Curves for this expression are shown in figure 3 
for values of L of 200, 600, 1,000, and 2,000 feet. At higher 
frequencies, these curves all approach a shape of l/fi 2 but 
differ in the frequency at which the flattening occurs. For 
increasing values of L, the curves flatten out at lower 
frequencies. Comparison of these curves with those of 
figures 1 and 2 arfd other measurements of the spectrum for 
atmospheric turbulence has indicated that the values of L 
for atmospheric turbulence are of the order of 1 ,000 feet. A 
value of 1,000 feet is used as a representative average value 
in the subsequent analysis. 


DERIVATION OF CONVERSION TECHNIQUES j 
BASIC METHOD 

In the present analysis, use will be made of some concepts 
and results in the theory of random processes. The theory 
of random processes is a recently developed branch of proba- 
bility theory and the aspects of the theory pertinent to the 
present applications are described in some detail in reference 
13. In particular, the relations between peak counts and 
spectra for a stationary Gaussian random process which 
will be used in the present study are derived therein. Some 
general aspects of random-process theory are also covered 
in references 20 and 21. 

In many recent studies of airplane behavior in rough air, 
atmospheric turbulence is generally considered a stationary 
Gaussian random process. The assumption that turbulence 
is a Gaussian random process appears warranted by the 
approximately Gaussian character of turbulent velocity 
fluctuations. The assumption of stationarity implies that 
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the statistical characteristics of the turbulence are invariant 
with space and time and also appears warranted for many 
purposes. For present purposes (in which the overall gust 
and load experiences of an airplane in operational flight are 
of concern), however, the process cannot be considered a 
simple stationary one, inasmuch as the turbulence charac- 
teristics of the atmosphere vary widely with weather con- 
ditions, particularly in regard to the intensity of the turbu- 
lence. 

In order to account for the variations in atmospheric turbu- 
lence with weather condition, it will be assumed that turbu- 
lence is only locally Gaussian and stationary; that is, its 
statistical characteristics are Gaussian and invariant in a 
given restricted region and for a short time but vary, particu- 
larly in intensity, from time to time and place to place. 
This assumption implies that the region or time is small 
relative to the entire flight path or flight duration but large 
enough for statistical equilibrium to be achieved. The 
Gaussian character for the turbulent velocity fluctuations 
for a given weather condition will be used as a building block 
to construct the operational gust history, which must 
cover many weather conditions and, in the overall, is not a 
stationary Gaussian process. On this basis, the overall tur- 
bulence experienced by an airplane in given operations is 
taken to consist of the summation for appropriate exposure 
times to a series of elemental stationary Gaussian processes. 

If the airplane response to turbulence is linear, as assumed 
in the present analysis, the response, such as the load 
history to each elemental turbulence process, is likewise a 
Gaussian process and the overall operational load history 
may in turn also be considered to consist of a summation of 
the loads for the various elemental Gaussian turbulence 
disturbances. This scheme serves to yield a reasonable 
approximation of the actual airplane load history and, 
further, has the particular advantage for present purposes of 
permitting the use of relations between peak counts and 
spectra derived for the stationary Gaussian case in reference 
13. 

The foregoing considerations form the basis for the present 
analysis. They will be applied in order to estimate the 
probability distribution of root-mean -square acceleration 
from operational data on peak accelerations. These esti- 
mates of the distribution of root-mean-square acceleration 
will then be used in order to obtain the associated probability 
distributions of root-mean-square gust velocity by taking 
into account the airplane dynamics. 

RELATIONS BETWEEN NUMBER OF PEAKS AND SPECTRA 

Simple Gaussian case. — The asymptotic relation between 
the average number of maximums per second exceeding a 
given value and the spectrum of a stationary Gaussian dis- 
turbance y(t ) has been derived in reference 13 and is for large 
values of y given by (see appendix) 



where 

N(y ) average number of maximums per second exceeding 
given values of y 
oo frequency, radians/sec 

$(w) power-spectral-density function of random disturb- 
ance y(t ) 



Equation (2) is the exact expression for the number of 
crossings per second with positive slope of given values of y 
but is an approximate expression for the number of maxi- 
mums above a given value of y. For present purposes, 
however, equation (2) appears to be an adequate approxima- 
tion of the number of peaks, as indicated in the appendix, and 
forms the basis for the present analysis. 

Examination of equation (2) indicates that the number of 
peaks per second above given values depends upon a 2 , which 
is the area under the spectrum, and upon the second moment 
of the spectrum about the origin. If it is assumed that the 
spectral shape of turbulence is invariant with weather con- 
ditions (as suggested by available measurements) and varies 
only in intensity or root-mean-square gust velocity, the out- 
put spectrum <U„(w) for acceleration a n (t) for a given airplane 
under given operating conditions (fixed speed, weight, 
altitude, etc.) is likewise invariant in shape. Under these 
conditions, the coefficient of the exponential term in equation 
(2) is a fixed constant and equation (2) may be written as 

N(a n )^Noe- a ^ 2 (3) 

where the constant No is given by 



The constant No gives the expected number of times per 
second that a n (t ) crosses the value zero with a positive slope. 
Because of the invariance of the shape of <fr a?i (w), the quantity 
N 0 is independent of the turbulence intensity. It has the 
dimensions of a frequency and can be considered a character- 
istic frequency of the airplane acceleration response to 
turbulence. 

If the logarithm of both sides of equation (3) is taken, the 
following relation is obtained : 

log iV(a„)=log iVo— ( 5 ) 

Equation (5) indicates that the log N(a n ) is a linear function 
of a 2 with slope equal to — l/2<r 2 . Thus, for a stationary 
Gaussian disturbance, the root-mean-square value a may be 
obtained simply from the slope of the line for the number of 
peaks when plotted on semilogarithmic paper as a function 
of a 2 . Because of the linearity of the relationship, this 
plotting arrangement is found useful. 
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A representative distribution of the number of peak accel- 
erations per second (both positive and negative peaks) ob- 
tained from operations of a transport airplane is shown as 
the solid curve of figure 4 as a function of a n 2 . The data are 
given, as is generally the case, for a threshold value of 0.3gf. 
The distribution appears concave upward and departs con- 
siderably from the straight lines that would be expected for a 
simple stationary Gaussian disturbance in this plot. It is 
thus clear that the overall distribution of peak accelerations 
cannot be adequately represented by the simple Gaussian 
case. A method for obtaining a more adequate representa- 
tion of the load history is explored in the following para- 
graphs. For this purpose, the idea of a combination or com- 
posite of Gaussian processes is used and equation (2), which 


expresses the relation between peak counts and the power 
spectrum for the simple Gaussian case, is extended to this 
composite condition. 

Composite Gaussian case. — If the overall operational load 
history (in terms of the average number of acceleration peaks 
per second exceeding given values) is considered to consist of 
various exposure times to different Gaussian gu§t disturb- 
ances, the average number of peaks per second exceeding 
given values is given by 

M(a n )—^iP iNi(a n ) (61 

i=l 

t 

where Pi is the proportion of total flight at the ith conditibn 
and Ni(a n ) is the number of peak accelerations per second 



Figure 4. — Illustrative example of graphical separation into Gaussian components of distribution of peak acceleration (operation 1). 
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exceeding given values of acceleration for the fth condition. 
In this form, equation (6) is general and permits accounting 
not only for variation in the spectrum of turbulence but also 
for variations in the airplane response characteristics. 

If equation (3) is substituted into equation (6), there is 
obtained 

M&=±| (N 0 )JP t e (7) 

i=l 

It will be assumed, as previously indicated, that the spectrum 
of atmospheric turbulence is invariant in shape but varies in 
intensity or in root-mean-square gust velocity. This assump- 
tion is suggested by the available measurements of turbulence 
spectra and appears reasonable for present purposes. Fur- 
ther, if the response characteristics for a particular airplane 
in operation are represented by a single-response function, 
such as obtained for average values of the airplane, and oper- 
ating parameters, such as the weight, airspeed, and density, 
equation (7) may be written for this (the discrete) case as 

Mfcj=ATo 2?^ (8) 

i= 1 

where the quantity No is, as in equation (3), fixed for a 
given airplane. 

Continuous case. — If the airplane is now assumed to 
encounter turbulence of all intensities (continuous variations 
in the root-mean-square gust velocity) but of the same spec- 
tral shape (fixed value of L for present applications), equa- 
tion (8) may be extended for the continuous case to yield 
the following relation : 

In this case, the sum of the terms of equation (8) becomes a 
continuous integral. The expression for the continuous case 
contains the same characteristic frequency No, which is inde- 
pendent of the turbulence intensity, the same exponential 
term as the earlier expressions, and, in addition, the function 
j(o a ^. This function is the counterpart of the ratios P x , P 2 , 

. . . P k in equation (8) and is the probability-density dis- 
tribution of c an which defines the proportion of time spent at 
various values of cr an ; the quantity f(o a ^do an gives the pro- 
portion of flight time spent at values of o an between o an and 

O' a -\-do a . 
a n 1 u n 

The relations between the number of peaks and the root- 
mean-square values of the random process given by equa- 
tion (8) for the discrete case and by equation (9) for the con- 
tinuous case are the basic relations in the present analysis. 
These two relations are applied in the following sections and 
methods are developed for estimating the distribution of the 
root-mean-square acceleration from operational gust-load 
data generally available in the form of peak counts M(a n ) . 
Then, the transformation of these distributions of root-mean- 
square acceleration to distribution of root-mean-square gust 
velocity is considered. 


DISTRIBUTION OF ROOT-MEAN-SQUARE ACCELERATION FOR DISCRETE 

CASE 

As previously indicated, the operational load history for a 
particular airplane in terms of the number of peaks .per 
second exceeding given values is given by equation (8). 
Each of the terms of the summation of equation (8) is given 
by 

( (10) 

which, as previously noted, yields a straight line of slope 
— 1 /2<r* 2 if log M(a n ) is plotted as a function of a n 2 . This con- 
dition suggests that the overall operational loads for a given 
operation when plotted in this specified form can be built up 
of straight-line or Gaussian components. A simple trial for 
a sample load history shown in figure 4 indicates that a good 
approximation to M(a n ) can be obtained with only a few 
components. (Note that the ordinate of figure 4 is 2 M(a n ) 
inasmuch as both acceleration peaks and acceleration mini- 
mums have been included in the operational data.) The 
values for the three components shown add up to give the 
short-dash curve which is seen to be a close approximation to 
the overall load history. The procedure devised for the deter- 
mination of these straight-line components consists of the 
following steps: First, line (1) is obtained by taking a tan- 
gent to the tail of the observed distribution; line (2) is then 
taken from the point on line (1) which underestimates the 
observed distribution by one-half and is drawn tangent to 
the upper part of the overall load-history curve. The third 
line, if required, is then obtained from line (2) in the same 
manner that line (2) was obtained from line (1). The sum 
of the values of the lines obtained in this manner will gener- 
ally yield a good approximation to the observed distribution, 
as will be seen subsequently. 

The procedure outlined in the preceding paragraph for 
obtaining the linear components of the observed distribution 
is somewhat arbitrary. Several alternative procedures were 
also considered and discarded. These procedures included 
the selection of the first component at the upper end of the 
curve and the use of specified combinations of slopes cor- 
responding to given values of root-mean -square gust veloci- 
ties (which might be considered averages for various weather 
conditions such as clear-air turbulence and thunderstorms). 
These alternative procedures appeared to offer no significant 
advantages and had the additional undesirable characteristic 
of yielding a poorer approximation at the larger values of 
acceleration which are of greatest interest. 

It will be recalled that the slopes of the lines in figure 4 are 
equal to — 1/2 <n 2 , —1/2 <r 2 2 , and — 1/2 <r 3 2 where o x , o 2 , and <r 3 
respectively, represent the root-mean-square acceleration for 
the three effective Gaussian components. The values for the 
case represented in figure 4 are o-j^O.430, o 2 — 0.247, and 
0-3=0.147. These values o u o 2 , and o 3 for the acceleration 
components are used subsequently to obtain the associated 
root-mean-square gust velocities for these Gaussian com- 
ponents. 

In order to estimate the values of the proportion of total 
flight time P { for the ith condition, it will be noted from 


471089 O -58 -2 
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equation (10) that for each component 


Pi 


Mj{a n ) 

N 0 e~ a nl 


(11) 


where Mi(a n ) is the number of accelerations per second 
exceeding a n for the fth condition or component. The num- 
erator of equation (11) may be obtained directly from the 
line for each component in figure 4. The denominator of 
equation (11), as would be expected, is the total number of 
peak accelerations per second exceeding a n that would be 
experienced if the airplane spent all its flight time at the 
indicated root-mean-square acceleration value. The value 
of the denominator is seen to depend upon the value of the 
constant No. 

The values of A T 0 and the denominator of equation (11) may 
be determined in a number of ways. They may be de- 
termined from their relation to the power spectrum as given 
in equations (3) and (4), but this procedure is not very 
practical for present purposes because the acceleration 
spectrum cannot be determined from the type of records 
available. The value of these quantities can also be de- 
termined in a number of ways from short sections of flight 
record in homogeneous rough air. \The most direct way is to 
make use of the property that No is equal to the average 
number of times per second that a n (t ) crosses the value zero 
with positive slope. Unfortunately, the iilm speed of the 
records available (2 to 8 feet per hour) was frequently too 
slow to permit such counts with the desired degree of 
reliability. 

Another method of determining the value of N 0 and the 
denominator of equation (11) from sample flight records 
involves the use of the property of equation (3) of defining 
the average number of peaks or crossings per second above 
given values of a n . If a n is set equal to ka in equation (3), 
the quantity Nika ) gives the average number of peaks per 
second that have a value of a n greater than ka where a is the 
value of the root-mean-square acceleration for the record. 
Thus, from equations (3) and (4) N(k<r), which is designated 
as N k for brevity, is given by 


1 

'J 

r* co — 

| w 2 $(u>)do) 

'o 

2 7T 

u 

<f>(w)do> 

'o 


yielding a root-mean-square acceleration above ka t is from 
equation (11) given by 


Pt 


M.ikot) 

' N k 


(13)' 


where M^kai) is the average number of peaks per second 
exceeding the value ka t for the linear component. This 
relation is used in the subsequent analysis of the operational 
loads data. In order to avoid the errors associated with the 
expression for the number of peaks at small values of k, a 
value of k of 2 was generally used. For the illustration pf 
figure 4, the values of 2M l {2a l ), 2M 2 (2<7 2 ), and 2M 3 (2o' 3 ) 
are indicated by ticks on the figure and are 4.3X10 -6 , 
3.2 X10 -4 , and 7.4 X10 -3 . The sum of the P’s will in any 
case be less than one, the remaining time being in either 
smooth air or very light rough air which does not contribute 
many peak values above the threshold value. 

DISTRIBUTION OF ROOT-MEAN-SQUARE ACCELERATION 
FOR CONTINUOUS CASE 

The representation of the load experience in the previous 
section in terms of a few discrete values of root-mean-square 
acceleration is a simplification, since atmospheric turbulence 
may actually be expected to cover a continuous variation 
in intensity. In this section, the problem of estimating 
the associated continuous distribution of root-mean -square 
acceleration from the overall peak counts will be considered. 

The determination of the actual probability distribution 
of root-mean-square acceleration f(a) from a given peak-load 
history requires the solution of the integral equation given 
by equation (9). Since it does not, in general, appear pos- 
sible to represent the load experience in terms of the number 
of peak accelerations above given values by a simple function , 
an effort was made instead to estimate /(a) directly in simple 
form, based on the results obtained in the application of the 
metnods for the discrete case. Equation (9) is then in- 
tegrated and the results obtained are preserted in chart 
form for use in comparison with operational data. 

Experience with the analysis of gust-load data, wherein the 
foregoing discrete-case methods have been used, has indicated 
that the distribution of a decreases rapidly with increasing 
values of a and might be approximated by simple exponen- 
tial distributions. On this basis, three exponential-type 
probability density functions were considered, namely, 


or 

Nk^Noa-^ ( 12 ) 

Inasmuch as the average number of peaks (or crossings) per 
second exceeding a value of a n ~ka can be estimated from 
the available sample records, equation (12) provides a 
useful method for determining the value of JVo and the 
denominator of equation (11). Care should be taken in 
such applications to insure that the record section used 
approximates a Gaussian disturbance. 

On this basis, the proportion of flight spent in rough air 


Case a: 




Case b: 

1 


y s (<r)=- ! 


a 2 

Case c: 

«')=sb* -V5/ - , 


( 14 ) 


where in each case cr^O. 
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For each case, a is a scale parameter, larger values of a 
representing more severe load histories. The three distribu- 
tions are shown in figure 5 for values of a= 1 and, in order, 
are seen to have increasingly larger areas at the higher values 
of o’. These functions, by the definition of a probability 
density function, all have unit area. 

Case a. — If it is assumed that f(o) is given by case a, 
equation (9) may be integ rated in closed form and yields the 
following result for M(a n )/N 0 : 


No 


(15) 


The function M(a n )jN 0 gives the proportion of peak acceler- 
ations exceeding given values of a n and is thus a cumulative 
probab ility distribution. Figure 6(a) shows the distribution 
M(a n )INo for a range of values of a x . Equation (15) yields 
a simple result for the proportion of peak accelerations per 
second exceeding given values of a n , in terms of the single- 
scale parameter of the distribution of The value of the 
par ameter u x is inversely proportional to the slope of the line 
for M{a n )jNo when plotted on semilogarithmic paper and 
also gives the value of o (ln below which 68 percent of the air- 
plane f light t ime is spent. The linear variation with a n of 
the log M(a n ) given by equation (15) is of considerable inter- 
est because gust-load flight-test data on peak accelerations 
under some conditions, particularly where a limited range 
of weather conditions is represented, tend to exhibit linear 
trends when plotted on semilogarithmic paper (log M(a n ) as 
a function of a n ). This condition, however, does not generally 
apply to operational data. 

Case b. — The distributions of M(a n ) observed from 
operational flights frequently exhibit a nonlinear trend and 



Figure 5. — Comparison of three assumed distributions. 


a concave upward form when the log M(a n ) is plotted 
against a n . This form implies that the distributions of 
f(o) for these cases might more adequately be represented 
by case b or case c. Substituting the expression for case b 
in equation (9) yields the following expression for the 
number of maximums exceeding a n : 

f "— « 2 ' 2 06 ) 

Jo 0,2 


A closed-form evaluation of the integral could not be found. 
Numerical evaluations of the integral, however, are made 
readily. If the substitution <r=sa 2 is made into equation 
(16) and both sides of the equation are divided by No, the 
following result is obtained: 


M(a n ) _ C 
No Jo 


a n 2 

6 42 * 


(17) 


Equatio n (17) was evaluated for various values of a« 2 /2a 2 2 
so that M(a n )jN 0 could be determined for various values of 
a 2 . The results obtained are shown plotted in figure 6(b). 
Operational gu^tjload data can be plotted conveniently in 
figure 6(b) in order to determine whether they adhere to the 
present distribution shape and, if so, to obtain the appropri- 
ate value of the scale parameter a 2 . (The plotting requires 
the determination of the value of No which may be estimated 
from flight records as described previously.) The points 
shown in the figure represent operational gust-load data 
which are discussed in a subsequent section. 

Case c. — For case c, equation (9) becomes 


M(a n ) 1 f” 

N o 2a 3 2 J o 


«n 2 _ Jg 

e ** 01 d* 


( 18 ) 


Equation (18) could not be evaluated in closed form but 
was evaluated numerically for various values of a 3 . The 
results obtained for a range of values of a 3 are shown in 
figure 6(c). Operational gust-load data can also be plotted 
conveniently in figure 6(c) in order to determine whether 
they adhere to the present distribution shape and, if so, to 
obtain the appropriate value of the parameter a *. 

In general, a good representation of the observed dis- 
tribution of peak values can be obtained by comparison of 
the data with the three charts of figure 6. 2 The choice of 
the appropriate case and the value of the scale parameter 
then specify the probability distribution of the root-mean- 
square acceleration in accordance with equations (14). As 
will be seen in the section entitled “Application to Gust- 
Load Statistical Data/’ these three distributions appear, in 
most cases, to be adequate for the representation of the 
types of data so far considered. However, other distribu- 
tion forms appear worth investigating in future studies. 


8 The charts of figures 6 (b) and 6 (c) contain extensions and revisions to the charts presented 
earlier in Technical Note 3362. 
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(a) Case a. 

Figure 6. — Cumulative probability of a peak -acceleration exceeding given values of a , 



Cumulative probability, 
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(b) Case b. 

Figure 6. — Continued, 






Cumulative probability, 
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DISTRIBUTION OF ROOT-MEAN-SQUARE GUST VELOCITY 

In the foregoing sections, methods have been given for 
converting the overall load histories from counts of peak 
accelerations into several discrete root-mean-square values 
of acceleration and their associated percentage exposure 
times to each level or into a continuous distribution of 
root-mean-square acceleration. The problem of the con- 
version of these results into the airplane gust experience is 
now considered. 

The root-mean-square gust acceleration in rough air is 
related to the turbulence spectrum and the airplane response 
characteristics by the following equation: 

= f $«,(“)&> (19) 

where 

$£/(«) power spectrum of vertical gust velocity 

T(u>) amplitude of airplane acceleration response to 

sinusoidal gusts of unit amplitude 
<i> are (w) power spectrum of airplane normal accelera- 

tion 

Since the root-mean-square acceleration obtained by the 
preceding evaluation is seen to depend on both the spectrum 
of turbulence and the airplane response characteristics, the 
root-mean-square acceleration is apparently not sufficient to 
fix the turbulence spectrum. Several procedures appear 
possible for the analysis of load measurements for the present 
purpose of deducing root-mean-square gust velocities. These 
include the actual measurement of acceleration spectra and 
the calculation of the airplane frequency-response function. 
The turbulence spectrum is then given by 

(«) 

° )= T^) (20) 

This procedure is extremely laborious and requires exten- 
sive calculations for the output spectrum the reli- 

ability of the results, in turn, depends upon the reliability of 
the calculated airplane transfer functions. At the present 
time, it is unlikely that the large amount of work involved in 
this approach is warranted even if the available records were 
in a form which would permit this type of analysis. Actually, 
the film speed used for the available records was far too slow 
to permit this type of analysis. 

As an alternative, some simplifications appear warranted. 
As a preliminary effort in this direction, it will be assumed 
that : 

(1) The airplane is rigid. 

(2) The airplane is free to move vertically only (not pitch.) 

(3) The spectrum of vertical gust velocity is given by 




L 1+312 2 Z 2 

7T (1 + fl 2 Z 2 ) 2 


These assumptions, although crude, should nevertheless 
provide some reasonable approximations of the gust histories. 
The possible magnitude of the errors introduced by these 
assumptions is considered in the section entitled “Reliability 
of Results.” 

For the foregoing conditions, a useful result obtained by 
Y. C. Fung (ref. 2) is that 4 


Z 2 = 


- <TU - - 


16V 2 
tcc\1+K)‘ 


I(K,s ) 


(21) 


where 

F 

c 

K 


mean-square acceleration 
mean chord 

airplane mass parameter, 


4 W 

girpSc 


s ratio of mean chord to scale of turbulence, c/L 

Since 1 in almost all cases of concern, equation (21) 
may be simplified to yield 


(Xu 


V' 


P VSm 1 1 (K, s) 
2 W 


(T a„ — A\Ct u 


( 22 ) 


where 


A? 


2W 


V" 


pVSm I(K,s ) 


and is an airplane gust-response factor defined in 

reference 2 and is shown in figure 7 as a function of K for 
various values of s—cfL. 

In this representation, the three-dimensional slope of the 
lift curve m has been used to replace 2ir in order to account 
for the overall three-dimensional aerodynamic effects. In 
this form, observed values of a 0n may be used directly with 
the values of the airplane parameters to determine <r v . 
This calculation requires the determination of the value of 
Ai. For this purpose, it would appear adequate to use 
average values of p, V, and W as well as an average value of 
K for a given operation or, if greater precision is required, 
for portions of an operation. Thus, the quantity A x may 
be treated as a constant for a given operation (or portion of 
an operation). 

It should be noted that the direct proportionality be- 
tween the root-mean-square gust velocity and acceleration 
as given by equation (22) is not restricted to the one-degree- 
of-freedom case but may be generalized to cover as many 
degrees of freedom as desirable. In general, 


<*a n=A<r v 

where, from equation (19), A is defined by 


(22a 


where 0 is a reduced frequency o>/V in radians per foot 
and L is the scale of turbulence. 


(22b 
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or for the gust spectrum assumed herein 

(22c) 

In these cases, the amplitude of the airplane frequency re- 
sponse function T(Q) will depend upon additional aerody- 
namic, geometric, and structural parameters of the airplane 
rather than upon those considered in the single-degree-of- 
freedom case. 

For the case of discrete values of o a n , equation (22) permits 
the direct evaluation of the associated values of root-mean- 
square gust velocity a v for the single-degree-of-freedom 
case. 

For the case of a continuous distribution of root-mean- 
square acceleration, the appropriate distribution of root- 
mean-square gust velocity is obtained from equation (22) or 
(22a) by the relation for a change of variables for probability 
distributions and in terms of /(o- 0w ) is given by 

f(*u) — AJiAiffu) {(r an =Ai<r v ) (23a) 

f (<T u) ~ Af (A<Tu) (<T an — A<T(j) (23b) 


Thus, for the three acceleration distributions considered, the 
distributions of root-mean-square gust velocity are given by 
the following equations: 


Case a: 


where 


Case b: 


where 


Case c: 


where 




, _ ff C/ 2 / 26 ! 2 


I <H 


k<?u)=f- e~°ul b 2 




b 2 —= 


'A i 


fz(<ru)= 2^2 « 


^ U3 

VA 


( 24 ) 
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DIRECT TRANSFER OF GUST LOADS FROM ONE AIRPLANE TO ANOTHER 

It may on occasion be desirable to transfer the gust loads 
directly from one airplane to another without recourse to 
first determining the probability distributions of root-mean- 
square acceleration and gust velocity. The following 
analysis indicates a simple direct manner of performing this 
transfer. 

Consider airplanes i and j. Assume that the gust loads 
have been measured on airplane i under fixed operating con- 
ditions and that it is desired to determine the load experience 
for airplane j under fixed operating conditions and with the 
same turbulence experience as that of airplane i. For 
airplane i, 


M(a)=(iVo), f“ da, (25) 

where the subscripts for acceleration have been omitted for 
brevity. Similarly, for airplane j 

f,(a,) e-’>W da, (26) 

If the turbulence experience is the same for both airplanes, 
then 


(27) 

Ai 

by virtue of equation (22) or (22a). Substituting equation 
(27) into equation (26) and making use of the relation for 
the change of variables for probability distributions yields 



The integral on the right-hand side is, from equation (25), 
equal to the quantity 1 


Thus, 



(29) 


(30) 


Equation (30) indicates that, in order to obtain the load 
histor}^ Mj(a n ), it is only necessary to take the value of 
measured distribution M { ( ) at the acceleration value 



teristic frequencies ( No)j/(No)i . Essentially, the operation 
involves only a change of the abscissa scale and a multiplica- 
tion of the ordinate value by the ratio of the characteristic 
frequencies. 


It should be noted that the values of A need not neces- 
sarily be restricted to the one-degree-of-freedom case con- 
sidered in the present conversion of the acceleration data tc 
gust velocities. In many cases, it may be desirable to’con- 
sider a more complete dynamic analysis in determining the 
quantities N 0 and A. 

6 

APPLICATION TO GUST-LOAD STATISTICAL DAT£ 

In this section, the methods developed in the section on 
derivation of conversion techniques are applied to available 
gust-load, statistical data in order to convert the statistical 
data in terms of peak counts of accelerations and gusts into 
distribution of root-mean-square gust velocity. The data 
to be considered are of two types: first, data obtained from 
a number of different transport operations, and, then, in ordei 
to obtain a description of the turbulence that is more flexible 
than the operational gust experience and that is applicable 
to arbitrary flight plans, available data on the variation oi 
turbulence with altitude and with weather condition. For 
the operational data, both the discrete and continuous 
methods are applied; whereas, only the continuous case is 
considered for the other data. 

OPERATIONAL GUST-LOAD HISTORIES 

Discrete case. — The graphical procedure outlined in an 
earlier section is applied to a series of eight operational loac 
histories obtained from NACA VGH records. The scope o: 
the data considered is summarized in table I. The datf 
were obtained from six different types of transport airplanes 
with operations 5, 6, and 7 representing the same airplant 
type flown by different operators. The pertinent airplam 
characteristics and operating conditions are given in tabk 
II. These data include estimates of the overall average 
flight altitude, average flight speed, and average weight 
The overall distributions of peak accelerations measured ir 
flight are given in table III. Because of the limited dati 
available, no effort was made to consider breakdowns of th( 
overall data by flight condition or altitude. Values of N, 
and N 2 , given in table IV, were estimated in accordant 
with equation (12) from sections of record obtained in con 
tinuous rough air by counting the number of peak accelera 
tions exceeding 2 a. These values contain some revisions o 
the values presented earlier in NACA Technical Note 336' 
but owing to the nature of the records available are stil 
considered crude estimates. The measured distributions o 
acceleration increments for the eight operations are showi 
in figures 8 (a) to 8 (h). (The illustration previously con 
sidered in figure 4 is operation 1 which is also shown ir 
fig. 8 (a).) The subdivision of the distribution into Gaussiar 
components in the manner previously described is also indi 
cated by the long-dash lines in each case. The values o 
a t and Pi for each operation obtained from the slopes of tin 
linear components and the values of 2M(2a t ) were deter 
mined in each case and are summarized in table IV. 


471089 0 - 58 -3 
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Figure 8. — Concluded. 
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The results of figure 8 and table IV indicate that the 
operational load histories can, in most cases, be reproduced 
by numbers of the order of 1 percent of the flight time at 
<r a? y=0.15 g and 0.05 percent of the flight time at a an =0.2>g. 
In some cases, three components were found desirable. 
Actually, the total flight time in rough air is considerably 
higher than the 1 percent indicated by these values, as 
previously mentioned, and should perhaps include as much 
as 10 percent at a lower value of a ajl of the order of 0.05(/. 
This flight time would yield the primary contribution to 
the number of peak accelerations at values of a n below 
0 . 3<7 but contributes only a few peak accelerations above 
the threshold value of 0.3 g used in the data evaluations 
considered herein. 

The conversion of the root-mean-square accelerations 
to root-mean-square gust velocities requires the deter- 
mination of the value of the quantity A x for each operation, 
as defined in equation (22), 


— pVSm 


V 


J(K,s) 


For this purpose, the airplane characteristics given in 
table II were used along with the average values of air 
density, airspeed, and airplane weight. In determining 
the gust response factor, average values of weight and air 
density were also used along with a value of L of 1,000 
feet. The values of A x obtained for the eight operations 
and the resultant values obtained for a v are given in table IV. 
It is of interest to note that a change in the value of L 
by a factor of 2 results in root-mean-square gust velocities 
that are increased by about 25 percent if L is doubled 
(£=2,000 feet) or decreased by about 25 percent if £ is 
halved (£=500 feet). 

Continuous case. — In order to determine the appropriate 
continuous distribution of root-mean -square acceleration, 
the tabulated distributions of table III were compared 
with curves of figure 6. Operations 1, 2, 3, and 8 appeared 
to follow case b best and are shown in figure 6 (b). (The 
measured distribution of peaks which included positive and 
negative peaks was divided by 2 for this comparison.) 
The appropriate values of a< z can be obtained from the 
figure by interpolation and are given in table IV. The 
distribution of <? a for these operations is, from, equations (14), 
given by 

/ 2 (Oa)'-=— - C-*al<h (<T ^ 0) 

Cb 2 


It is of interest to note that all three of the two-engine 
aircraft in low-altitude operations, operations 1 , 2, and 3, 
are best represented by this case. 

The remaining four operations had a slower rate of decrease 
in M{a n )JN 0 with increasing a n ; this condition suggested 
that case c might be more appropriate. The data for these 
operations are shown plotted in figure 6 (c) and are in fair 
agreement with the curves for case c. Again, one-half the 
measured distributions are plotted to be comparable with 


the calculated curves. The acceleration histories for these 
operations are thus described by the probability distribution 
for case c (eqs. (14)) 

where the appropriate value of a 3 determined By inter- 
polation from figure 6 (c) for each operation is given in 
table IV. 

In some of the foregoing cases, the acceleration histories 
are not too well represented by the analysis. Better repre- 
sentations, if desirable, can be obtained, however, by con- 
sidering combinations of the three distributions or more 
complicated distribution forms. 

The conversion of these distributions to distributions of 
root-mean-square gust velocities only requires the deter- 
mination of the values of b defined in equations (24). The 
appropriate values of b (b 2 or b 3 ) were determined by use of 
the values of A x for each operation and are given in table IV. 
These values apply to the associated distributions which are 
given by equations (24). The probability distributions 
obtained for each of the operations are shown in figure 9 (a). 
Because of the changes in the values of No and several other 
minor changes, these results also differ in some cases from 
those presented in NACA Technical Note 3362. The 
cumulative probability distributions, which are obtained by 
integrating the probability distributions, are given in figure 
9 (b) and define the proportion of total flight time spent in 
turbulence exceeding given values of g v . 

The description of the gust experience in this form is 
directly applicable to load calculations for other airplanes in 
similar operations by reversing the procedures used in obtain- 
ing these results. However, direct application of these 
results would apply only to similar operations. In order to 
obtain results that are more flexible and applicable to arbi- 
trary flight plans, it would be desirable to determine the 
variations in these distributions with altitude, weather con- 
dition, and, perhaps, geography. Preliminary results ob- 
tained in this respect are described in the next sections. 

VARIATION WITH ALTITUDE 

Conversion of data.-— In order to arrive at some rough 

A 

estimates of the variation oif(a v ) with altitude, use was made 
of the summary of gust statistics given in reference 22. 
Figure 6 of reference 22 presents estimates of the average 
gust experience at various altitudes that are representative 
of contemporary types of transport operations. In order to 
estimate the associated distributions of root-mean-square 
gust velocity, these results which are in terms of the fre- 
quencies of derived gust velocities were, for convenience, 
first converted to acceleration peaks by using the conven- 
tional gust-response equation as given, for example, in refer- 
ence 12, and the characteristics of a representative transport 
airplane. The charts of figure 6 could then be used to 
estimate the appropriate distribution form and the scale 
value a. 
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(a) Probability density. (b) Cumulative probability. 

Figure 9. — Probability distribution of the root-mean-square gust velocity for eight operations. 


The airplane and operating parameters used for this conver- 
sion to accelerations are given in table Y. The acceleration 
histories as defined by the cumulative probability distribu- 
tion of peak accelerations M(a n )jNo obtained on this basis are 
shown in figure 10 (a) for each 10,000-foot altitude bracket 
with the exception of the lowest altitude bracket. Also 
shown in the figure are curves for M(a„)/No corresponding to 
the probability distribution of root-mean-square acceleration 
/ 3 (<r a ) as given in figure 6 (c). It will be noted that, except 
for the lowest altitude bracket, the acceleration histories can 
be reasonably well approximated by members of this family 
of distributions. The appropriate values of a 3 obtained by 
interpolation on the figure for these altitudes are given in 
table V. For the lowest altitude bracket, 0 to 10,000 feet, a 
better representation was needed. A few trials indicated 
that a very good approximation for the data for the lowest 
altitude bracket could be obtained by the following two-term 
probability distribution : 

/( % )=(0.99) (31) 


The distribution of peak accelerations M(a n )INo correspond- 
ing to this distribution for a dn is shown in figure 10 (b) as the 
short-dash curve and is seen to be in excellent agreement with 
the data points for this altitude bracket. The contribution 


of each of the two terms of the right-hand side of equation (3F 
to the total curve is indicated by the two long-dash curves. 

By applying the values of A x for each altitude bracket ir 
table V and the foregoing results in equation (23a), the follow' 
ing distributions /(o-f/) were obtained for the various altitudt 
brackets : 

For the altitude bracket of 0 to 10,000 feet, 

/(<r„)= 0.99 jig <TV‘-«+ 0.01 g-ij e-'v» M (32 


for both altitude brackets of 10,000 to 20,000 feet and 20,00( 
to 30,000 feet, 


K'u> 


2(0.32) ; 




(33 


and for both altitude brackets of 30,000 to 40,000 feet anc 
40,000 to 50,000 feet, 

M ' 




2(0.29) : 


0 Vv 0 * 29 


(34 


In both equations (33) and (34) average values of b 3 ar< 
given for the two 10,000-foot intervals covered. These dis 
tributions are shown in figure 11 (a). The -cumulative dis 
tributions are shown in figure 11 (b) and give the proba 
bility of exceeding given values of <r v . 





Cumulative probability, 
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Acceleration, g units 


(a) 10,000 to 50,000 feet. 

Figure 10. — Probability of a peak acceleration exceeding a given value of a„ for a representative transport in flight at various altitude brackets. 




Cumulative probability, 
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(b) 0 to 10,000 feet. 
Figure 10. — Concluded. 
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(a) Probability density. 

Figure 11. — Distribution of <tu wi 


Perhaps the most important points to be noted from 
figure 11 (b) are the relatively large amount of time spent 
in essentially smooth air at the higher altitudes (<r J/ <C2 feet 
per second 93 to 96 percent of the time) and the relatively 
large amount of time spent in light to severe turbulence at 
the lowest altitude bracket (<^^>2 feet per second 25 percent 
of the time). The time spent above 5 feet per second for 
the lowest altitude bracket is roughly five to ten times as 
great as that for the higher altitude brackets. 

It should be remembered that these results are in terms of 
true gust velocity. If equivalent gust velocities which are 
more directly related to the airplane response are used, the 
decrease in turbulence intensity with altitude would, of 
course, be even more pronounced than is indicated herein. 

Method of application.— The method of application of 
these results to the calculation of response histories for 
arbitrary flight is based on the following relation: 

_ t In 2 

MW=JV.riW« (35) 

which is obtained by substituting equations (22a) and (23b) 
into equation (9). The procedure involves the division of 
the operational flight plan into homogeneous portions or 
segments in regard to flight altitude and operating con- 
ditions such as airspeed and airplane weight. The appropri- 
ate distribution of the root-mean-square gust velocity is 



selected for each flight segment from figure 1 1 (a) and equa- 
tion (35) is then evaluated for each segment. Actual numer- 
ical calculations are facilitated by the use of the charts of 
figure 6. The sequence of steps involved in such applica- 
tion is as follows: 

(1) The operational flight plan is divided into homogene- 
ous segments in regard to flight altitude (10,000-foot alti- 
tude brackets) and operating conditions such as airspeed 
and weight. 

(2) The appropriate distribution of j{<r v ) is selected for 
each altitude bracket from figure 11 (a). 

(3) The values of A are determined for each significant 
segment of the flight plan in accordance with equation (22) 
or (22c). 

(4) In order to obtain the associated distributions of accel- 

A 

eration f(a a ), each of the distributions of is transformed 
by the relation 

/ (°'a„) = i f(*u) 

A 

where 

(Tu—ffajA 

(5) The values of iV 0 are most easily determined from 
flight records, if available, by the methods already indi- 
cated. For new designs N 0 must be estimated analytically. 
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From equations (4) and (20) 


and for thunderstorms, 


iVo= 


2t 



o> 2 $ f/ (w)T 2 (co)c?co 

4> f/ (aj)7 T2 (co)rfco 


I** 


The reliable estimation of No from this relation is a diffi- 
cult problem inasmuch as the value depends heavily on the 
reliability of both the gust spectrum and the airplane fre- 
quency response function at the higher frequencies. In 
order to obtain reliable estimates of this quantity, it appears 
necessary to account for the effects of spanwise variations in 
turbulence which act strongly to attenuate the airplane 
response at these higher frequencies. The effects of the 
spanwise variations of turbulence on the airplane response 
have been studied in references 6 and 23 and will have to 
be considered further in order to establish how best to 
account for these effects in estimating the value of N 0 by 
analytical means. 

(6) The distributions of a a>i and the values of N 0 are then 
used in equation (9) to derive the number of peak accelera- 
tions per second or per mile for each condition. These cal- 
culations are facilitated by the use of charts such as given 
in figure (6). 

(7) The results obtained in step 6 are then weighted in 
accordance with the flight distance in each condition or 
segment and then summed for all conditions in order to 
obtain the overall acceleration history. 

VARIATION WITH WEATHER CONDITION 

In the preceding discussion, the variation in the gust ex- 
perience with flight altitude was considered. Another 
breakdown of the gust experience which may be useful in 
some problems is the variation in gust experience with 
weather conditions. Figure 12 (a) shows estimates of the 

A 

variations in /(oy) that have been obtained from data for 
several types of turbulent weather conditions (ref. 24) . The 
curve for clear-air turbulence was based on data obtained in 
flight through clear and turbulent air at the lower altitudes. 
The curve for cumulus clouds was based on data obtained 
in flight under moderate convective conditions such as rep- 
resented by bulging cumulus clouds. Finally, the curve for 
thunderstorms was based on data obtained in flight in the 
immediate vicinity of or within severe thunderstorms. In- 
asmuch as the data of reference 24 showed a linear variation 
of the logarithmic gust frequency with gust intensity, the 
probability distributions of a v is given by case a (fig. 6 (a)). 
The curves of figure 12 (a) are thus given by the equation 
for case a: 

e -^ 2 /^ 2 ( au ^ 0 ) 


The values obtained for b x are as follows : 
For clear-air turbulence, 

&i = 3.15 

for cumulus clouds, 

&i=6.28 


6i= 10.05 

A simple linear measure of the relative intensity of the tur 
bulence for these conditions may be obtained by comparing 
the values of &i. It will be noted that the values of b t for tin 
cumulus and thunderstorm conditions are roughly two anc 
three times that for the clear-air condition. The cumulative 
probability distributions for the three weather conditions are 
shown in figure 12 (b). 

It has been estimated that contemporary transport ,opera 
tions spend about 10 percent of their flight time in this clear- 
air-turbulence condition, 1 percent in the cumulus condition 
and perhaps 0.05 percent in thunderstorms. These results 
may therefore be applied in evaluating the effects on the over- 
all load experience of operational procedures which woulc 
tend to modify this weather experience. For example, the 
introduction of airborne radar for weather avoidance may b( 
expected to reduce the exposure time to the severe turbulence 
conditions of thunderstorms. Also, high rates of climb anc 
descent through the lower and more turbulent altitude layers 
may cause a drastic reduction in the 10 percent of the ffighl 
time attributed to clear-air-turbulence conditions and there- 
by cause a marked reduction in the overall gust experience 

RELIABILITY OF RESULTS 

The foregoing analysis has indicated that, under reasonable 
assumptions, available gust statistics in the form of counts ol 
acceleration and gust peaks may be converted into a forir 
suitable for applications in spectral calculations of airplane 
responses. The appropriate form for this purpose is the prob- 
ability distribution of the root-mean-square gust velocity 
Estimates of this distribution are derived from data for £ 
number of transport operations, and the variation of this 
distribution with altitude and weather condition is also given 
In the derivation of these results, a number of assumptions 
were made and these assumptions may affect the reliability ol 
the results obtained and their method of application in re- 
sponse calculations. The reliability of the results and the 
method of applying these results are discussed in this section 

The principal errors in the present analysis stem from the 
following sources: 

(1) The assumption of the form of the turbulence spectrum 
and the, assumed value of 1,000 feet for the scale of the 
turbulence 

(2) The assumption of average values of weight, altitude 
and airspeeds in the reduction of the data 

(3) The restriction of the airplane dynamics to the one- 
degree-of-freedom case in the conversion of the results from 
acceleration to gust velocity 

(4) The assumption that the gusts are uniform across the 
span 

(5) The statistical sampling errors that arise from the 
limited data considered. 

The magnitude of the errors resulting from each of these 
sources is considered in turn. 
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(a) Probability density. (b) Cumulative probability. 

Figure 12. — Distribution of root-mean-square gust velocity with weather condition. 


Turbulence spectra.— The assumed spectral shape given 
by equation (1) appears to be a reasonable representation of 
the spectra measurements reviewed in figures 1 and 2 for 
the region of gust wavelengths covered (10 feet to perhaps 

10.000 feet). At the shorter gust wavelengths, the spectrum 
probably approaches zero at a faster rate than that assumed. 
However, the contribution of these higher frequencies to the 
total power is very small and appears negligible. 

At the very low frequencies (say, gust wavelengths greater 
than 3,000 feet), the spectrum is not yet adequately defined 
although the available measurements suggest a flattening of 
the spectrum in this region and values of L of the order of 

1.000 feet. It may be expected that the value of L varies 
somewhat with weather condition and, as previously indi- 
cated, such variations would affect the reliability of the de- 
rived root-mean -square gust velocities. This assumption 
thus appears to be a possible source of error in the distribu- 
tions of root-mean-square gust velocities derived herein. 

The consequences of such possible errors in the scale of 
turbulence and in the root-mean-square gust velocities on the 
reliability of load and other response calculations will vary 
depending upon the airplane response characteristics and 
have to be evaluated in each case. For many airplanes, the 
errors in the calculated loads may be largely compensating 
inasmuch as the errors in the values of root-mean-square 
gust velocity will be balanced by the error due to the use of 
the erroneous value of L (L= 1,000 feet) in the calculation of 
the quantity A. This condition can be verified in a concrete 
fashion for the one-degree-of-freedom case by some trial cal- 
culations by use of the curves for the gust response factor of 


figure 7 and consideration of possible variations in the value 
of the scale of turbulence. For example, representative 
values of the mass parameter and the mean chord for the 
airplanes used to collect the gust data are 100 and 10 feet, 
respectively. For these values and the scale of turbulence 

c 

of 1,000 feet, the value of is 0.01 and the value of the 


gust response factor 


ig 0 465< 


If the true value of 


the scale of turbulence is 2,000 feet, the value of J-JEH: 

about 0.36. The use of the higher value of the gust response 
factor would, from equation (22), yield a root-mean-square 
gust velocity that is about 29 percent too low. However, the 
determination of the value of A 1 for the airplane for which the 
loads are being calculated would, in turn, be high by about 
the same amount for many cases. For example, for the same 
mean chord but with a mass parameter of 200, the value 
obtained for the gust response factor in figure 7 is about 0.61 
for L— 1,000 feet, as compared with the value of 0.48 for the 
scale of turbulence of 2,000 feet. Thus, the use of a value of 
L= 1,000 feet in the load calculation yields a compensating 
error of about 27 percent. Thus, at least for limited varia- 
tions.in the airplane parameters, the net errors in the load cal- 
culation due to variation in the scale of turbulence from the 
assumed value of 1,000 feet would be- small. 

Average values of airplane parameters— The effects of 
assuming average conditions of weight, altitude, aiid speed 
may also be expected to introduce some errors in the derived 
root-mean-square gust velocities given in figure 9 for the 
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various operations. The errors resulting from the assump- 
tions of average weight and altitude may he expected to he 
small because these errors should largely average out. The 
effects of speed variations are in a somewhat different cate- 
gory inasmuch as efforts are normally made to reduce speed 
in rough air. For the operations considered herein (fig. 9), 
it appears that the reductions in speed from normal operating 
speeds were generally small and negligible at the lighter levels 
of turbulence. At the more severe levels of turbulence, the 
airspeeds were, on the average, somewhat lower than normal 
operating speeds, but the reductions in most cases were 
small, about 5 to 10 percent. Thus, the distribution of a v 
at the higher levels may be biased to this extent. 

In deriving the results on the variation of the distributions 
of root-mean-square gust velocity with altitude and weather 
condition, figures (11) and (12), the assumption of an average 
speed was not necessary inasmuch as the original distribu- 
tions of derived gust of reference 8 were largely based on 
data in which the actual speeds were taken into account. 
However, average conditions of altitude and weight were 
assumed both in the derivation of the basic data and in the 
present reductions. The errors due to these sources should 
be largely self-balancing and therefore small. 

Airplane dynamics. — In the derivation of the root-mean- 
square gust velocities, the airplane was assumed to be rigid 
and restrained in pitch although free to move vertically. 
These assumptions are admittedly rough approximations 
and their effects on the derived gust data require additional 
study. Extensions of the present results to include the effects 
of these two additional degrees of freedom appear possible 
although they would involve a considerably larger number of 
airplane and operating parameters. It should, however, be 
possible in the meantime to make some rough overall cor- 
rections to the root-mean-square gust velocities for some of 
these effects. For example,, available information suggests 
that the neglect of the dynamic structural-response effects 
on the center-of-gravity accelerations might be expected to 
have resulted in roughly a 10-percent over estimation of the 
root-mean-square gust velocities for most of the airplanes 
considered in the present study. The effects of the airplane 
pitching motions on the root-mean-square gust velocities for 
the present airplanes are also generally considered small. 
Recent studies suggest that the pitching motions for the 
present airplanes would tend to decrease the gust accelera- 
tions and thus tend to lead to some underestimation of the 
root-mean-square gust velocities. A rough estimate of 
about 10 percent appears reasonable for this effect for the 
airplanes considered. Thus, it is suggested that the overall 
effects of pitch and flexibility might largely cancel each 
other. However, additional study of these effects is needed. 

Effect of spanwise variations in turbulence. — In the 
present analysis, no consideration has been given to the 
the effects of gust averaging introduced by the finite span. 
If the ratio of the span to the scale of turbulence is large, 
these effects, as has been indicated in reference 6, may be 
appreciable. For example, for a span of 150 feet and a 
value of L= 1,000 feet, the results of reference 6 suggest that 
the root-mean-square value of c v may be underestimated 
in the present analysis by perhaps as much as 5 to 10 percent. 
For smaller values of the ratio of the span to the scale of 


turbulence, the magnitude of this discrepancy is consider- 
ably smaller. If desirable, the incorporation of these effects 
on an average basis would also appear straightforward. 

Statistical reliability of results. — Since the present results 
are based on relatively small operational samples (of the 
order of 1,000 hours), the statistical reliability of the desired 
distributions of root-mean-square value is dependent upon 
the statistical reliability of the acceleration data. Past 
experience has indicated that these distributions are reliable 
for samples of this size at the lower levels of acceleration 
(0.3<7 to 0.5gr) but have poor statistical reliability at the 
higher acceleration values. As a consequence, it may be 
expected that the derived probabilities for the higher root- 
mean-square gust velocities are only rough estimates and 
should be used only as a guide. More reliable information 
regarding the higher root-mean-square gust velocities 
requires more extensive flight data although it may be 
possible to supplement the present results by use of available 
NACA V~G records. The extension of the present appli- 
cations to include such other data is, however, beyond the 
scope of the present report. 

Applicability of results. — In view of the foregoing con- 
siderations, the distribution of root-mean-square gust veloc- 
ity given herein may be considered a reasonable first-ordei 
estimate of the characteristics of atmospheric turbulence 
that are essentially independent of the characteristics of 
the airplanes involved. Thus, these gust spectra and 
root-mean-square gust-velocity distributions can be reason- 
ably applied in gust-load calculations in which the effects ol 
pitching motions and flexibility are included in the deter- 
mination of the airplane transfer functions. 

A final word of caution is appropriate in regard to tht 
limitations of the present results, particularly figure 11. 
for response calculations. Inasmuch as the basic data or 
which the present results are based were largely obtainec 
from conventional American transport-type operations 
they are representative of such operations in regard to sucl 
factors as geography, turbulence avoidance procedures, ant 
terrain clearance. As a consequence, the present result; 
may require modification for operations that differ in i 
significant fashion from those considered herein. 

CONCLUDING REMARKS 

The foregoing analysis has served to demonstrate thai 
the gust statistics may, under reasonable assumptions, b< 
converted into a form appropriate for spectral-type calcula 
tions. The significant and fundamental quantity, for thii 
purpose, appears to be the probability distribution of th< 
root-mean-square gust velocity. The results obtained ii 
defining the variations of this function with type of operation 
flight altitude, and weather condition provide at least i 
starting basis for their application to response calculation! 
in arbitrary operations. These results should serve tc 
supplement the discrete-gust techniques in current use ant 
be particularly appropriate in problems requiring a mor< 
detailed accounting for the airplane dynamics than is 
possible by discrete-gust techniques. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 5, 1956. 



APPENDIX 

RELATIONS BETWEEN PEAK COUNTS AND SPECTRA FOR A GAUSSIAN RANDOM PROCESS 


In the present analysis, use is made of the relations be- 
tween the number of maximums per second and the spectra 
for a Gaussian random process. These relations are derived 
in reference 13 and are summarized herein in order to per- 
mit the examination of the reliability of the approximate 
expression used in the body of the report. 

NUMBER OP MAXIMUMS 

The probability p m that a Gaussian random process y(t) 
will have a maximum intensity ranging from y x to y x + dy x 
in the time interval t to t-\-dt is given in reference 13 as 

p m =dyidl e— m„», , / 2 |m| + 

*i *<2 kJ e ~ n H 1+ erf dSbO] (A1) 

where the error function is 


erf Z x = 


_2_ rz x 

V^Jo 


e~ z2 dZ 


and the coefficients M i} can be expressed in terms of the value 
at zero of the autocorrelation function and its deriva- 


tives as 

M n = — * o'Vo <4) 

Af„=(* o") 2 

Mj,%~ Mo" 

|Mh-^Wo (4) -(W 2 ] „ 


(A2) 


where 


+( T )= l jZ 2 l ff T y(t)y(t+r)dt 


and \p o" and i^ 0 (4) are, respectively, the values at 0 of the 
second and fourth derivative of the autocorrelation function 
i p(r). The values at zero of the autocorrelation function 
and its derivatives are in turn related to the power-spectral- 
density function <f>(co) by the following relations: 



Equation (Al) is rather un wieldly, but two simple results 
of interest may be obtained from it. The first result— the 
expected number of maximums per second N p — is obtained 
26 


by integrating over y x from — » to °° and over t for a time 
of 1 second to obtain t 

N,=~ W-o 

zir 



(A4) 


The second result is an asymptotic expression for large y x 
for the number of maximums per second N(y x ) exceeding 
given values of y X} which can be obtained by integrating an 
asymptotic approximation to p m from ?/i to ; the result, 
given in reference 13, is 





It is of interest to note that the right-hand side of equation 
(A5) is also the exact expression for the number of crossings 
per second with positive slope of given values of y x . 

Equation (A5) is the basic relation used in the present 
analysis for the number of peaks per second exceeding given 
values of y x . Since it is an approximation for this purpose, 
the magnitude of the errors introduced in the present analy- 
sis by its use is of interest and is considered in the remainder 
of the appendix. 


RELIABILITY OF THE APPROXIMATION 

Past experience has indicated that in gust-load applica- 
tions, the approximation given by equation (A5) is in most 
cases good for values of y x /<x^> 2. For values of 2 /i/<r< 2 , 
equation (A5) tends to underestimate the number of peak 
loads to some extent. The magnitude of these errors does 
not, however, appear to be large and, as will be indicated, 
has only a very small effect on the reliability of the present 
analysis. 

The magnitude of the error introduced by the approxima- 
tion for small values of y x /a may be indicated by considering 
the ratio N(0) to the total number of peaks N p . From 
equations (A4) and (A5), this ratio is given by 


m o) J 

r- 

'o 

2 4>(a>) dco 

N P , 

(J 

f a) 4 <t(a>) (foy 
'o / 
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For a low-pass filter, which appears to describe roughly the 
response of a relatively rigid airplane to turbulence, equation 
(A6) reduces to 


vV(0)_y5 ^ 
N p ~ 3 


(A7) 


The quantity N p , which gives all the maximums, includes 
some maximums at negative values of y x . For the low-pass- 
filter case, the results of reference 13 indicate that about 15 
percent of all maximums are at negative values of y t . Thus, 
the approximation of equation (A5) appears to be roughly 
10 percent low for the low-pass-filter case at y x /<T—0. For 
increasing values of y x j a, this error decreases rapidly and is 
less than 3 percent at yi/<x= 0.5. 

For moderately flexible airplanes of the type considered in 
the present study, the degree of underestimation of the 
asymptotic formula is somewhat larger than for the band- 
pass case. In this case, calculations indicate that equation 
(A5) appears to underestimate the number of peaks by 
about 30 percent at y x ja— 0, 10 to 15 percent at y 1 /<r=l, 
and 2 to 3 percent at yi/<r= 2. The effect of these errors is, 
however, considerably mitigated in the present applications 
for the following reason. The total number of peak accel- 
erations exceeding given values is seen from equation (9) 
to depend in principle upon the whole distribution of 
root-mean-square values. Actually, for the particular 
exponential-type functions for f(a) conside red in the present 
analysis, the principal contributions to M(a n ) arise from 
values of a n ja an that range from about 1.0 to 4.0. Thus, 
the asymptotic expression is being applied principally over 
the region where the underestimation is only a few percent. 
Therefore, the errors introduced in the present analysis by 
the use of the asymptotic formula (eq. (A5)) can be con- 
sidered negligible. 
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TABLE I.— SCOPE OF OPERATIONS 


■ 

Operation 

Route 

Flight 

hours 

Flight 

miles 

i 

Northern transcontinental .. 

834 

186, 000 

2 

Rocky Mountains — North and South... 

331 

55,500 

3 

Southern transcontinental 

676 

130,000 

•1 

90 percent east of Mississippi River 

771 

202,000 

5 

New York to Europe — New York to South 
America. . 

1,079 

284,000 

6 

San Francisco to Honolulu 

1,953 

486,000 

7 

Northern transcontinental 

876 

235,500 

8 

Southern transcontinental. ... _ ... .... 

706 

193,500 


TABLE II— AIRPLANE AND OPERATIONAL CHARACTERISTICS 


Operation 

Gross 

weight, W, 
lb 

Average 
weight, W, 
lb 

Wing 
area, S, 
sq ft 

Average 
chord, c, 
ft 

Average 
flight 
altitude, ft 

Average air 
density, p, 
slugs/cu ft 

Airspeed, 

V, 

ft/sec 

m, 

per radian 

Mass 

parameter, 

K 

£=1,000; 

V™ 


39,900 

33, 915 

864 

10.1 

5,000 

0. 002049 

327 

5.0 

75.00 

0.411 


25,200 

21,420 

987 

10.4 

8,000 

. 001869 

246 

4. 92 

44. 14 

.318 


40, 500 

34, 425 

817 

9.7 

5,000 

.002049 

281 

5.03 

83.83 

.426 


107, 000 

90, 950 

1,650 

14.7 

10,000 

. 001756 

384 

4.93 

84.55 

.385 


142,500 

121, 125 

1,720 

12.9 

12, 500 

. 001622 

386 

5.12 

133.08 

.567 


142,500 

121, 125 

1,720 

12.9 

12,500 

.001622 

366 

5.12 

133.08 

.567 


142,500 

121, 125 

1, 720 

12.9 

12,500 

.001622 

394 

5.12 

133. 08 

.567 


89,900 

76, 415 

1,463 

13.7 

12,100 

.001653 

402 

4. 95 

91.20 

.495 









TABLE III— NUMBER OF ACCELERATION PEAKS 
EXCEEDING GIVEN VALUES 



Cumulative frequency for various operations 1 


1 

2 

3 

4 

. . 

5 

6 

7 

8 

0.3 -- 


19,483 
4, 632 

5, 593 
1,350 
365 

1.888 

659 

612 

909 

1,287 

232 

0.4 


152 

132 

178 

0.5 

Bin 

1,288 

iKf?S 

40 

31 

40 


0.6 

377 

370 

104 


13 


14 

18 

0.7 

124 

100 

31 

17 

5 

4 

5 

6 

0.8. . ... 

47 

35 

16 

6 

2 

2 


1 

0.9 

26 

9 

5 

2 

1 

1 



1.0 

7 

3 

2 


1 



1.1 

4 


2 


1 




1.2 

2 








1.3 

1 
















Total flight hours.. 

834 

331.1 

676.5 

770.8 

1,078.5 

1, 953. 4 

875.5 

706.5 


1 Number Includes both positive and negative peaks. 

















A REEVALUATION OF DATA ON ATMOSPHERIC TURBULENCE FOR APPLICATION IN SPECTRAL CALCULATIONS 2 { 


TABLE IV— RESULTS FOR VARIOUS OPERATIONS 

(a) Summary of acceleration experience 



Acceleration experience for operation — 


1 

2 

3 

4 

5 

6 

7 

8 8 

<Ta 

0. 430 

0.287 

0. 323 

0.278 

0.349 

0.364 

0. 255 

0.226 

<y a 

0.247 

0. 181 

0. 181 

0. 147 

0. 178 

0. 171 

0. 132 

0. 128 


0.147 

4.3X10-6 
3.2X10-* 
7.4X10-3 
1.0 
0. 270 
1.59X10-5 
1. 18X10-3 
2. 74X10-2 
0.030 





0. 104 



<r ° n 3 

2Mi(2<n) 

2M 2 (2<7 2 ) 

2Afo(2cro) 

No 

2 V 2 
Pi 
P 2 
Ps 
02 

03 

1. 35X10-* 
6.2X10-3 

1.35X10-5 
8 . 5X10-* 

1.9X10-5 
8 . 2X1 O'* 

9. 5X10-1 
6 . 6X10-5 

4.6X10-2 
3. 7X10-5 
3.0X10-* 

9. 5X10-6 
5. 4X10-* 

2. 85X10-5 
1. 15X10-3 

1.0 
0.270 
4. 99X10-* 
2. 29X10-2 

1.0 
0.270 
5.00X10-5 
3. 15X10-3 

0.5 
0. 135 
1. 41X10-* 
6.07X10-3 

0.5 
0. 135 
7. 04X10-6 
4. 89X10-* 

0.5 
0. 135 
3. 41X10-6 
2. 74X10-* 
2. 22X10-3 

0.5 

0.135 

7.04X10-5 

4.00X10-3 

0.5 
0. 135 
2. 11X10-* 
8. 52X10-3 

0.034 

0. 026 




0.022 

0.041 

0.036 

0. 034 

0.037 









(b) Summary of gust experience for £=1,000 ft 



Gust experience for operation — 

1 

2 

3 

4 

5 

6 

7 

8 

Ai 
a i'i 
«v 2 

62 

63 

0.01755 
24. 501 
14. 074 
8. 376 
1.709 

0. 0166 
17.289 
10.904 

0. 0146 
22. 063 
12.363 

0.0146 
19. 041 
10. 068 

0.0129 
27. 033 
13. 788 

0. 0122 
29. 787 
13.993 
8.511 

0. 0132 
19. 347 
10. 015 

0. 0156 
14. 515 
8. 221 

2.048 

1.781 




1.410 

.339 

0.317 

0.308 

6. 322 






TABLE V— AIRPLANE AND OPERATIONAL CHARACTERISTICS USED FOR DISTRIBUTION OF a v WITH ALTITUDE 


[ IF— 76,415 lb; ,8=1,463 sq ft;~=13.7 ft; Vi=401.5 ft/sec; ?«= 4.95 per radian; Vo =0.5 per sec] 


Altitude, ft 

Air density, 
p , slugs/cu ft 

Mass 

parameter, 

K 

h (K,,) 

» 7 r 

Ai 

d2 

«3 

&2 

63 

O tn 10V103 

0. 002049 

73.4 

0. 451 

0. 01758 

0. 026 


1.479 



.050 


2.844 


10 to 20 

. 001496 

100.5 

.524 

. 01491 


0.040 


0.327 

20 to 30 

. 001065 

141.2 

.583 

.01181 


.035 


.313 

30 to 40 

. 000736 

204.5 

.659 

.00907 


.0275 


.289 

40 to 60 

.000458 

328.4 

.741 

.00638 


.0225 


.282 
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